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Abstract 

In a model proposed in 2012, all flavor mixing has a single source and is governed by a single 
“master matrix.” This model was shown to give several predictions for quark and lepton masses and 
mixing angles and for mixing angles within SU{'5) multiplets that are observable in proton decay. 
Here it is shown that the same master matrix controls the flavor-changing processes mediated by 
a singlet scalar that exists in the model, giving predictions for r —>■ /ry, r —>■ ey, and /i —>■ ey. . 


1 Introduction 

In a 2012 paper [T], we proposed a model in which all flavor changing effects, including CKM mixing 
[1] and MNS mixing [3], are controlled by a single “master matrix.” In that paper, the model was 
shown to give several predictions for neutrino properties, including the Dirac neutrino CP phase, as 
well as post-dictions for quantities that are still not precisely known, such as the atmospheric and 
solar neutrino mixing angles, and ms/rrid- In a subsequent paper [2] it was pointed out that the same 
model predicts all the mixing angles that come into gauge-boson-mediated proton decay, thus giving 
further tests. 

In this paper, we show that the same model gives predictions for flavor-changing effects produced by 
the exchange of a Standard-Model-singlet scalar field that exists in the model. In principle, therefore, 
certain parameters of the model could be measurable in three independent ways: by precise determina¬ 
tion of neutrino and quark properties, by proton decay branching ratios, and through flavor-changing 
decays such as r —)• /ry, r ey, and /U —)• ey. 

The model is based on two assumptions: (1) that SU{5) symmetry relates quarks and leptons, 
and (2) that all flavor violation comes from mixing between three chiral fermion families that we 
shall denote lOj + 5i, i = 1, 2, 3, and N vector-like fermion multiplets that we shall denote 5^ -|- 5^, 
A = 1 ,..., N. (N can be as small as 2.) In particular, it is assumed that the Yukawa terms that involve 
only the three chiral families are flavor-diagonal due to an abelian family symmetry group, which we 
shall call Gf- ■ The vector-like fermion multiplets, on the other hand, do not transform under Kp, 
and as a consequence their mixing with the three chiral fermion families produces flavor violation 
among those families. More specifically, the mixing responsible for flavor violation is between the 5* 
and the 5^^ multiplets. This means that the model is of the “lopsided” type [5]. Because the origin 


of flavor changing is in the 5 sector, it shows up more strongly in left-handed leptons (which are in 5 
multiplets) than in left-handed quarks (which are in 10 multiplets). This gives a simple and elegant 
explanation of the fact that the MNS angles are much larger than the CKM angles, as is the basic 
idea of so-called “lopsided models” [5]. 

The effect of the mixing of Sj and the 5^ shows up in the effective low-energy theory of the 
known quarks and leptons as a 3 x 3 non-diagonal “master matrix,” which we call A, that appears 
in their mass matrices. By field and parameter re-definitions this master matrix can be brought to a 
simple triangular form (which we call ^a), which contains only one complex and two real parameters, 
whose values can be completely determined from CKM mixing. This allows predictions for all other 
flavor-changing effects. 

In section 2, we will review the model and show how it leads to predictions for flavor changing in 
the lepton sector and in proton decay. In section 3, we will analyze the flavor-changing effects that 
arise from the exchange of a Standard-Model-singlet scalar that exists in the model. 


2 Review of the model 


We shall now review the model and its predictions for masses and mixing matrices. More details can 
be found at mm- The Yukawa terms of the model are 


J^Yuk = Y*(10a0*)(5/^) + 2/*(10A)(5^) 

+ Y,{10ilO^){45H) +U^0A){45 ^h) 
+ iXi/MR)(Si5i){5H){5H) 

+ YXb{5a5b){1h) +y'A^i5A5^){rHi), 


( 1 ) 


where the subscript H denotes Higgs multiplets. Repeated indices {i, A, or B) are summed over. The 
first two lines contain typical Yukawa terms that give realistic quark and lepton masses. The third 
line is the effective Weinberg dim-5 operator that gives mass to the neutrinos in either the type-I or 
type-II see-saw mechanisms. All the terms in the first three lines involve only the multiplets 10 * -|- 5 * 
and are therefore flavor-diagonal. The fourth line of Eq. (1) contains the Yukawa terms that give 
mass to the vector-like fermions and mix 5 * with 5 ^. These masses, coming from S'C/(5)-singlet Higgs 
fields, can be much larger than the weak scale, and indeed can even be of order the GUT scale. The 
only assumption that is required to fit the CKM and MNS mixing angles is that the masses generated 
by these two terms are of the same order, which we shall refer to as the “heavy scale” M*. This scale 
must be large enough to explain why these new fermions have not been observed. The fermions that 
do not get mass of order M*, which consist of the 10 * and three linear combinations of 5 * and 5a, are 
the known quarks and leptons, which we will call the “light fermions”. 

There are many abelian family symmetries that could enforce the flavor-diagonal form of the terms 
in the first three lines of Eq. ( 1 ). A simple example (though not the simplest) is that given in [T], 
namely Gp = Ki x K 2 x K 3 , where (for a given i equal to 1, 2, or 3) Ki is a Z 2 symmetry under which 
10 *, 5 *, and I'^i are odd and all other fields even. Note that the vacuum expectation values of 1'^ 
spontaneously break the abelian family symmetries; so that the last term in Eq. (1), which mixes the 
5 * and 5 ^, does not respect the family symmetries and can give flavor violation. It is important for 


the predictivity of the model that the last two terms in Eq. (1) involve only 5?7(5)-singlet Higgs helds, 
as otherwise the “master matrix” would be different for quarks and leptons. This can be ensured by 
another abelian symmetry that prevents the SU{5) adjoint Higgs field from coupling in these terms 

[I]- 

The Yukawa terms in the first three lines of Eq. (1) give rise to the following diagonal mass 
matrices niu, nrid, m£, and rrii,- 


'^i Uj, di {iTid^ij dj, £i £j, Vi {Tni£)j uj ( 2 ) 

These are not the mass matrices of the known fermions, because we have not yet taken into account 
the mixing of the 5 multiplets, of which N linear combinations are “heavy” and 3 are “light”. A 
block-diagonalization to separate the heavy and light 5 fermion states is needed in order to find the 
effective mass matrices for the light fermions. 

Let us call [Y^^(l//)] = Mab and [UAii^Hi)] = ^Ai- Let us first examine the down-type quarks. 
These have a (3-|-A^) x (3 + A^) mass matrix of the form 



^ rrid 

V ^ 


0 

M 


V J 


( 3 ) 


The block diagonalization is carried out by a bi-unitary transformation of the (3 + N) x (3 + N) mass 
matrix: 


^ nid 0 

^ A M 



/ Gt 




^ rrid 

V ^ 


0 
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A B 
\C D 


^Ur 


( Md 

V 0 



( 4 ) 


Here the elements of G are small and Ul is approximately diagonal, because the elements of are 
very small compared to those of M and A. One can give exact expressions for the matrices A, B, C, 
D, and G, which will be useful in section 3. Defining T = M“^A, one can write 


A = [/ + rtr]-V2^ 

B = = Art = rt[/-hrrt]-V2 = 

G = -T[I + T^T]-^/^ = -TA = -[I + TT^-^/^T=-DT (5) 

D = [/ + rrt]-V2 

G = 

Since the elements of A and M are of the same order, the elements of T are of 0(1), and the matrices 
A, B, G, D have off-diagonal elements of 0(1). By simply multiplying out Eq. (4) one sees that the 
effective mass matrix of the three “light” down-type quarks, namely Md, is given by 


Md = nidA. 


( 6 ) 












The reason that nid gets multiplied on the right by the matrix A is that the matrix rud originally 
appears in a term 10i{md)ij5j, and the matrix A represents the mixing of the 5 multiplets. The mass 
matrix m£ of the charged leptons appears in a term 5i(m£)ijl0j, as can be seen from Eq. (2), and 
so gets multiplied on the left by . Therefore, the effective mass matrix of the three light charged 
leptons is Mg = A^nii. The up quark mass matrix appears in a term which involves 

no 5 multiplets, and so does not get multiplied by any factors of A. Therefore, Mu = mu- Finally, 
the mass matrix of the neutrinos, which comes from the dim-5 Weinberg effective operator, appears 
in a term 5i{mi^)ij5j, and so gets multiplied on both the right by A and the left by . Hence, we 
have altogether 


Mu = mu, Md = mdA, Mi = AAm^, My = A^myA. (7) 

We thus see that all flavor violation is controlled by A. Moreover, the matrix A can be brought to a 
simple form in the following way. By multiplying A on the right by a unitary matrix, the elements 
below the main diagonal of A can be made zero. Then by rescaling the rows by multiplying from the 
left by a complex diagonal matrix, the diagonal elements of A can be set to 1. That is, A can be 
written 


A = V A^U, 


( 8 ) 


where 2? is a complex diagonal matrix, 7/ is a unitary matrix, and A^ is a matrix of the form 




I ^ ^ 
0 1 

^ 0 0 



(9) 


where a, b, and c are real. It is easily seen that the matrix U can be absorbed into redefined right- 
handed down quarks and the left-handed lepton doublets. Similarly, the phases in T> can be absorbed 
into redefined fields. The diagonal real matrix \'D\ can be absorbed into redefinitions of the original 
diagonal mass matrices as follows; rfld = md\T>\, m£ = m£|D|, ffly = my\'D\‘^, and = mu- Thus, 
after these redefinitions, the mass matrices of the three light families take the new form 


Mu = mu, Md = mdA^, M^ = A'^m^, My = A'^myAA- ( 10 ) 


It is easy to see that to a very good approximation the elements of the diagonal matrix nld are just 
the eigenvalues of Md, i-e- the physical masses of the d, s, and b quarks. Therefore, in the basis of 
Eq. (10), the mass matrices of the up quarks and down quarks look as follows 


Mu 


^ mu 0 0 \ 

0 me 0 ; 

0 0 mt J 


/ 

md 

0 

0 \ 


/ 

1 

b 



/ 

md 

bmd 

ce md \ 


0 

ms 

0 



0 

1 

a 

= 


0 

ms 

ams 

V 

0 

0 

mb ^ 


1 

0 

0 

1 J 


1 

0 

0 

mb J 


( 11 ) 


One sees immediately that 



( 12 ) 


IT/ 1 am.s _ 

1 mb 


iT^ 1 ^ bm^ 

1 Vus\ — rns ~ 


jr . ~ ce'-^mg _ 
rrih 



0 = Skm- 


Similarly, the elements of the diagonal matrix mi in Eq. (10) are to a very good approximation the 
masses of the e, /r, and r. In the basis of Eq. (10), therefore, the charged lepton mass matrix Mf has 
the form 


/ 1 0 0 


Mf ^ 


ce 


id 


1 

a 


/ me 0 0 

0 m^ 0 

0 0 itIt 


me 

bme 

JO 


m. 


y ce me am^ 


0 

0 

mr 


(13) 


This is not diagonal, but the rotations required to diagonalize it are very small for left-handed charged 
leptons (namely, of 0(m^/m^), 0(m^me/m^), 0{m‘l/m‘^)). Thus, to a very good approximation, in 
the basis where the charged lepton mass matrix is diagonal, the effective neutrino mass matrix 
has the form (from Eqs. (9), (10) and (12)) 



1 

0 

0 ■ 


■ qe^f^ 

0 

0 ‘ 


■ 1 

ms 1 T/ 1 
us\ 
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0 

pe®“ 

0 


0 

1 

^\yco\ 



^.\ycb\ 
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0 

0 

1 


0 

0 

1 


We have scaled out an overall mass scale and parametrized the diagonal matrix rflv as diag(qe*^, pe*", 1). 
There are nine neutrino observables: three masses, three MNS mixing angles, the Dirac CP-violating 
phase, and two Majorana CP-violating phases. These are determined by five model parameters, p, q, 
a, 13, and Therefore there are four predictions, which we may take to be {My)ee (which comes 
into neutrino-less double beta decay), and the three CP-violating phases. In [1], it is found that the 
model’s best-fit values are 5Dirac = I-IStt radians, and {M^)ee = 0.002 eV. 

But the model actually is considerably more predictive than counting parameters suggests, due 
to the fact that the expressions for observables in terms of model parameters are very non-linear. It 
is found that for good fits, certain quantities that have already been measured (such as, 6atm, dgoi, 
and ms/md) must be in a restricted part of their present experimental range. Eor example, a value of 
the atmospheric angle smaller than vr/d is preferred by the fits, and a value of mg/md less than the 
median value of 20 is somewhat preferred. (See [T] for details.) As shown in [2] the model also makes 
non-trivial predictions for branching ratios in proton decay, which we will not review here. 


3 Flavor Changing from Singlet Scalar Exchange 

In this section we consider the effects of the scalar field \h that couples to the vector-like fermions to 
produces the N x N mass matrix Mab = '^ab{^h)- We will henceforth call this singlet Higgs field 
H = (H) -|- H. The exchange of the H will mediate flavor-changing processes. Eor these effects to be 
observable in practice, we must assume that the scale M*, which characterizes the mass and vacuum 












expectation value of 17, is not too much larger than the weak scale. We will assume that it is of order 
1 to several TeV. 

Let us look first at the Yukawa couplings of 17 to the down-type quarks. In the same notation of 
Eq. (3), the Yukawa couplings of 17 to the down-type quarks is given by 


(<^(10)) 11(5')) 

^ 0 0 \ 



y 0 M/(17) ) 

V Ut) ! 


(15) 


When one block-diagonalizes to separate the light and heavy fermion stats, this Yukawa matrix is 
transformed by the unitary matrices Ul and Ur as in Eq. (4): 

^0 0 \ / 

V 0 M/(17) ) [ 

So the effective Yukawa coupling of 17 to the three light down-type quarks d, s, and b, is given by 
di{G^MC)ijd'^ (17/(17)). Remarkably, this Yukawa matrix, which we will call can be written simply 
in terms of the master matrix A. Using Eq. (5), one gets 


I Gt 
-G I 


=ul 


/ 0 0 
^ 0 M/(17) 


( A B 
^ G D 

^Ur 


1 


/ G^MG G^MD 

y MG MD 


• (16) 


Y^(17) ^ G^MG ^ 


{-mdT^D‘^M-^)M{-TA) 

nidT^D'^TA 

nidT^I + TT^)-^TA 

mdT^T{i -h rir)-M 

md{A~‘^ — I)A^ = md{A — AA"^ A). 


(17) 


In going from line 3 to line 4, we have used the fact that (I -|- TT'^T)~^T = T{I + T^T)~^, as can 
easily be seen by expanding out the expressions is parentheses as power series. In the last line, we 
have used the fact that A is hermitian. Let us rewrite this expression in terms of the triangular matrix 
A^, since that is the matrix whose elements are known. Using Eq, (8) we have 


Yrf(17) = md{A - AA^A) 

= mdiVAAU - {VAaU){U^A^^V*){VAaU)] (18) 

= mdVAA[I - A^^\V\‘^Aa]U. 

The factor U on the right will be absorbed by the re-definition of the right-handed down-quark fields 
that was discussed after Eq. (9). Doing this re-definition, and using the fact that nidVAA = = 

Md, the Yukawa coupling matrix takes the form 


Yd{n)=Md[I-Al\VfAA]. (19) 

The mass matrix Md is diagonalized by a bi-unitary transformation to give V^MdVR = = 

diag(mrf, rus, mb). From Eq. (11). One sees that the matrix Vr is the CKM matrix, while the matrix 



Vr differs from the identity matrix by terms of order 0{ml/ml), 0{mdms/ml), and which 

can be neglected. It is clear then that in the physical basis of the down quarks 

^ - A^^IV^Aa]. (20) 

{it) 

Obviously, only the second term in the brackets leads to flavor changing. Let us parametrize the 
unknown matrix D as diag((5, e, ()). The flavor-changing Yukawa coupling matrix of the fl to the 
physical down-type quarks is of the form di{Yf'^)ijdji}, where 
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( 21 ) 


where 


Ad. = Asd = \5\H, 

Ad, = A^, = |<5|2ce*^ 

A,, = Al = \e\^a + \6\Hce^^. 


( 22 ) 


Note that the flavor-changing (i.e. off-diagonal) elements of YJ'^ depend only on two unknown 
combinations of parameters: |(Ip/(n) and |ep/(n). Note also that Ad. and A.d are real in the 
physical basis of the quarks, so that the e parameter of the system does not put constraints 

on flavor changing coming from the singlet scalar exchange. 

The charged-lepton sector is identical except for a left-right transposition. So writing the flavor¬ 
changing Yukawa coupling matrix of the to the physical charged leptons as if {Yf'^)iji~(t, one 
hnds 
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(23) 


where 


Ae^, = A^e = \5W 

Aer = A%e = \5?Ce^f 

A^r = A%^ = \efa + \5fbce'-^. 


( 24 ) 



The flavor-changing Yukawa couplings come into the processes ^27 through-two loop di¬ 

agrams, as shown in [6]. The specific diagrams that dominate in this model have the vector-like 
fermions running around the loop that gives an effective Q-photon-photon coupling. The resulting 
branching ratios for the flavor-changing lepton decays can be expressed in terms of the quantities 
given in Eq. (24) as follows [7j: 


BR{h £ 27 ) = 24 




(25) 


One prediction is that 


BR{t ey) ^ 

Aer 
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C 

to 

1 

cr 

to 

Vub 

BR{pL ey) 

Ae/i 
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KrnsJ 

Vus 


(26) 


If one assumes that the expression for A^r in Eq. (24) is dominated by the |(5p term, then one would 
also have the prediction 


BR{t —/ry) = \c\‘^BR{fi —)• ey) = 16 • BR{fj. —)• ey). (27) 

Given the present limit [8] that BR{fi ey) < 5.7 x 10“^^, this would gives a prediction that 
BR{t —^uy) < 10“^^. This is well below even what is expected to be observable at a super-c-r 
factory [9]. On the other hand, the branching ratio for this decay can be much larger if in Eq. 
(23) is dominated by the |ep term. As we will show below, there is an approximate theoretical bound 
that |ep < 1/2. This would give 


/1 Tp V \ 4 

i?i?(r ^ ^y) < 1.5 X 10-9 . (28) 

The flavor-changing processes involving quarks do not get large enough contributions from the 

exchange of the singlet scalar to stand out from Standard Model contributions. Eor instance, 

2 

the coefficient of {sd){sd) operators is found from from Eqs. (23-24) to be of order < 

( \ 4 ^ 

I , where we have used an upper bound on that is derived below. (From the 

first line of Eq. (32) one finds that [dp must be less than (1 -|- ~ 1/17-) 

Let us now consider the parameters d, e, (. While the matrix R = diag((5, e, C) is not known a 
priori, it is nevertheless possible to derive strict upper bounds on the parameters |(5|, |e|, and |C| from 
the properties of the master matrix A. From the fact that A = {I + and that A = RA/^JA, 

one has that 


A At = pAaA^p* = {i + rtr)-i 

(pAaA^p*)-^ -1 = Ttr. 


( 29 ) 


Computing the matrix on the left side of the above equation, one obtains 














I = T^T. (30) 


1/|5|2 -b/{S*e) (a6-ce*®)/((5*C) 

—b/{6e*) (l + 6^)/|ep —(a + a6^ — 6ce*®)/(e*C) 

(a6 — ce“*®)/((5C*) — (a + a6^ — 5ce“*®)/(eC*) (1 + + |a6 — ce*®P)/|CP 

For any matrix T, there is an inequality that must be satisfied by the elements of T^T. namely 

|(rtr),,f <(TtT),,(TtT),y, V z,j. (31) 

This is obvious if we write Tij = where = 1,2,3, are three complex vectors. Then the 

inequality is just seen to be the statement that • t(j)| < |t(j)|. From this inequality with 

{i,j) = (1,2), (1,3), and (2,3), respectively, one obtains after a little algebra 

(l + 62)|5|2 + |e|2 < 

(l + a2 + |a6-ce*'^|2)|,5|2 + |C|2 < l + a^, (32) 

(1 + + |a6 — ce*^P)|ep + (1 + 6^)|CP < 1 + 6^ + c^. 

using the values of a, b, c and 6 given in Eq. (12), the third equation of Eq. (32) gives an upper bound 
on |ep of approximately 1/2, as used in Eq. (28). 

In this paper, we have assumed that the scale M* of (12) is in the low TeV range, as otherwise 
the flavor-changing effects from exchanges of 12 would be hopelessly small. But then one must run 
the Yukawa couplings Y^g and shown in the last line of Eq. (1) from the GUT scale down to 
the scale M*. If these ran differently for the leptons and quarks, it would make the matrices A and 
M in Eq. (3) different for quarks and leptons, and thus also make the master matrix A different for 
quarks and leptons. That could destroy the predictivity of the model. If one considered only gluon 
loops in the running there is no problem, as the effect would be to increase A and M by the same 
factor for quarks relative to leptons. This factor would cancel in the ratio T = M~^A, and therefore 
also in yl = [/ -|- However, the gluon loops could do the following: they could increase the 

Yukawa couplings and for quarks to such an extent that the effect of these Yukawas on their 
own running could be much more significant for quarks than for leptons. That would make the forms 
of the matrices A and M — and therefore the form oi A — different for quarks and leptons. 

There are two ways to avoid this problem. One is that all the Yukawas Y^^ and remained small 
for the whole range from Mqjjt to M*. This has a drawback, however. If these Yukawa couplings 
Y^g are small compared to 1, then the VEV (12) would have to be large compared to a TeV to make 
the vector-like fermions in 5' -|- 5 heavy enough not to be seen. That would suppress flavor-changing 
effects from 12 exchange. 

A cleaner way to avoid the problem is to assume the following two conditions; (a) The Yukawa 
couplings that generate the mass matrix A are small compared to 1, and the VEV of the Higgs 
fields l( 4 j correspondingly large compared to a TeV. (That would have the additional advantage of 
making flavor changing from the exchange of these scalars negligible.) (b) The Yukawa coupling matrix 
^AB proportional to the identity matrix, which could be the result of a flavor symmetry that acted 
on the vector-like families. Then even if gluon loops drove Y^^ to be of order 1 at low scales, that 
would not affect the form of Y^^. 



Another theoretical issue raised by M* being near the weak scale is that the spontaneous breaking 
of the family symmetry group Gp would cause cosmological domain walls. This breaking is done by 
the VEVs To avoid overdosing the universe, these domain walls would have to be “inflated 

away”. One simple possibility is that Gj? is actually broken at a scale much higher than M* but only 
induces a VEV for Ipi that is of order M*. For example, consider the terms C{a) = +ipil;a, 

where the scalar field a and fermion bilinear are odd under a Z 2 and M is of order the GUT scale. 
Let the fermion bilinear get a condensate (V'V') = where A ~ which is many orders of 

magnitude bigger than M*. The Z 2 will be broken at the scale A, whereas {a) = A^/M^ ~ M*. 

4 Conclusions 

The model of flavor symmetry and flavor violation proposed in [1] has the virtue that it is (a) concep¬ 
tually simple, (b) explains some of the qualitative features of the quark and lepton spectrum {e.g. the 
MNS angles being much larger than the CKM angles), and (c) is highly predictive. As such, it can 
provide a kind of “benchmark” for seeing how large various kinds of flavor-changing processes might 
be expected to be. 

The model is of the “lopsided” type [5], which tends to give relatively large flavor-changing effects. 
In models with symmetric mass matrices, which are very common in the literature, off-diagonal Yukawa 
couplings Yij are typically proportional to yfmpmjlv. In lopsided models, however, Yij and Yji are 
very different in magnitude from each other, being proportional to nii/v and mjjv. This is the reason 
for the name “lopsided”, and why the flavor-changing effects tend to be relatively large. 

It is likely, then, that the flavor-changing Yukawa couplings given in Eqs. (21-24) (with the bounds 
in Eq. (32)) are typical of what would expect for a new scalar field. We see that if the scale of new 
physics M* is of order 1 TeV, there is good hope of eventually seeing the processes r —^ 7 , r ^ eq, 
and /i — 7> eq. One also sees from this model, that observing such processes can confirm or rule out 
specific models of the origin of flavor and flavor violation. 
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